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ABSTRACT 

Top quarks can be produced from the electroweak Wtb coupling, leading to a 
single top or antitop in the final state. We examine single top quark production at the 
Fermilab Tevatron and calculate the cross section as a function of top quark mass, 
for each of three separate production processes. We give details of the interferences 
between Feynman diagrams in W-gluon fusion and then study the effects of an 
anomalous (V+A) term in the Wtb coupling. Using predictions for the next Tevatron 
run, we estimate the experimental sensitivity to the CKM matrix element V t b as a 
function of the strength of the anomalous coupling. 

1. Introduction 

The top quark has recently been discovered by the D0 1 and CDF 2 collaborations at 
the Fermilab Tevatron. Measurements of its mass are 199 + 30 GeV and 176± 13 GeV 
from D0 and CDF in the lepton+jets decay modes and 145+32 GeV in the dilepton decay 
mode from D0 3 . Because the top quark is so heavy, electroweak production of single top 
quarks from a Wtb vertex becomes competitive with strong production of tt pairs from a 
gtt vertex. The resummed next-to-leading order cross section for tt production 4 is 
2.26^ i7 pb at 4s =1.8TeV when the top quark mass m, = 200 GeV. A recent 
calculation to 0(a 2 a s ) for single top production 5 gives the cross section for the principal 
single top process (for t and t combined) as 1.7 pb, which is 75% of the total it rate. 
Therefore, the possibility exists for an experimental measurement of single top production 
at the Tevatron with the current -100 pb 1 data set. 

We have examined the processes for single top production, and calculated the 
corresponding cross sections, at the current Tevatron energy of 4s =1.8 TeV. In addition, 
we have started to prepare for the large numbers of single top events predicted for the next 
Tevatron run at 4s = 2.0 TeV, when there should be twenty times as much data as exists 
now. We have examined the effects of an anomalous right-handed (V+A) coupling at the 
Wtb vertex from which single top quarks are produced. The single top quark system is the 
only place to study this coupling, and to measure the CKM matrix element V tb . 

There is an extensive literature on single top physics at hadron colliders 5 ' 6 , although 
few papers use the current Tevatron energy with modern parton distributions. How to 
study the Wtb coupling and what can be learnt have also featured in the literature 7 . 
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2. Single Top Quark Cross Section 



2.1 Single Top Processes 

We have calculated the tree level cross sections for the following processes: 



1. pp^tb + X 

2. pp^tq + X 

3. pp-^tW~ + X 



where X represents additional final state particles other than gluons. In process 1. the final 
state t and b are produced from an extremely off-shell s-channel W boson. Process 2. 
occurs via a t-channel W boson. Process three has several types of Feynman diagrams. We 
have included the following subprocesses in our calculations: 



where q is a valence or sea u or d quark. Subprocess 2.2 g'g — > tqb is known as W-gluon 
fusion. Although 1.2 and 2.2 look superficially the same, in fact they each have two 
Feynman diagrams which form separate gauge invariant sets, and the diagrams are 
calculated in pairs as higher order corrections to the 2— >2 processes with tb and tq in the 
final state. Feynman diagrams for these processes are shown in Fig. 1. 

In our calculations, we have omitted subprocesses with strange or charm sea quarks 
in the initial state for simplicity. These would boost the total cross section by -2-3%. We 
have included diagrams with intermediate state photons or Z bosons where appropriate (3.2 
and 3.3) but the contributions are not significant. For subprocesses 3.2 and 3.3, we have 
omitted diagrams with strong tt production, where the i has decayed to a Wb . We have 
included diagrams with off -diagonal CKM matrix elements. 

The subprocesses we have included in our calculations comprise all the significant 
ones with two or three vertices, except those with a gluon in the final state, which are 
significant but have not yet been fully treated. Subprocesses with an extra quark in the final 
state, for instance bg^tqq' and qb^tWq, although they have several Feynman 
diagrams, do not contribute more than -1-2% to the total cross section. This also applies to 
the subprocess bb — > tWb , despite its large number of Feynman diagrams, including one 
with electroweak tt production. 

2.2 Calculation Details 

We have calculated the production cross section for each of the subprocesses 
mentioned above. We used the computer program CompHEP 8 to do the tree level symbolic 
calculations and to generate optimized FORTRAN code for the squared matrix elements. We 
used the BASES 9 package to integrate over all phase space using parton distribution 
functions, and a CompHEP-BASES interface 10 to generate the event kinematics with 



1.1 q'q^>tb 
2.1 q'b^tq 
3.1 bg^tW 



1.2 q'g^tbq 
2.2 q'g^tqb 
3.2 qq^tWb 



3.3 gg->tWb 
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ja) 1.1 ud^>tb 1.2 ug-itbd 




jc) 3.1 bg -» W 



Fig.l Representative Feynman diagrams for the significant subprocesses included in 

our calculation for single top quark production at the Tevatron: (a) W boson s- 

channel production pp —> tb + X; (b) W boson t-channel production 
pp -» tq + X and (c) pp -> + X . 

smoothing of singular variables. For parton distributions, we have used the CTEQ3M 11 and 
MRS(A') 12 fits, which are representative next-to-leading-order sets in the MS 
renormalization scheme. 

The following Standard Model parameters have been used in our calculations: m z = 
91.19 GeV, m b = 5.0 GeV, a = 1/128, sin 2 tf w = 0.225 and CKM matrix elements V lld = 
0.975 and V, h = 0.999. All results have been obtained in the unitarity gauge and the 
't Hooft-Feynman gauge, as a check of the calculations. Differences between calculations 
in the two gauges are less than 0.1%. 

We have chosen to use mf as the QCD evolution parameter or scale (£. A typical x 
value is 0.1, where x is the fraction of the proton or antiproton momentum carried by each 
initial state parton. At a scale £ 2 = (180 GeV) 2 , the value of a, is 0.104. 

2.3 Combining Cross Sections 

Care must be taken when combining some single top subprocesses to avoid double 
counting. Subprocesses 2.1 (cfb—ttq) and 2.2 (q' g — > tqb : W-gluon fusion) have a 
considerable region of overlap when the g^bb is on-shell, so that the b quark which 
fuses with the W boson is the same as the initial state b sea quark in subprocess 2.1. A 
technique has been developed 13 which takes care of this double counting whilst being 
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consistent with the definition of the b sea quark in the parton distributions. It involves 
summing the cross sections calculated for each of the two subprocesses, and then 
subtracting the gluon splitting rate convoluted with subprocess 2.1. We have adopted this 
technique here. 

The subtracted term is not a small correction. For example, using CTEQ3M and m, = 
180 GeV, the 2— >2 subprocess q'b^tq gives a cross section of 0.750 pb before 
subtraction, W-gluon fusion q'g^tqb gives 0.293 pb and the correction is -0.536 pb, 
giving a total rate for pp — > tq + X of 0.507 pb. The subtraction term forms just over half 
of the raw total rate. Alternatively, if the 2— >2 subprocess and the subtraction term are just 
ignored, then the remaining cross section from only the W-gluon fusion subprocess is less 
than 60% of the total pp — > tq + X rate. 



2.4 Cross Section Versus Top Quark Mass 

Our results are shown in Fig. 2 as a function of the top quark mass m, , at ~Js = 
1.8 TeV. The solid curves are from the calculations with CTEQ3M and the dot-dash curves 
are from MRS (A'). With m, = 180 GeV, we find the total cross section for top to be 
0.834 pb (CTEQ3M), with contributions of 0.247 pb from subprocess 1.1, 0.032 pb from 
1.2, 0.215 pb from 2.1 after correction, 0.293 pb from 2.2 (VF-gluon fusion), 0.042 pb 
from 3.1, and 0.005 from subprocesses 3.2 and 3.3 combined. 




3. A Closer Look at W-Gluon Fusion 

We have studied the contributions to the production rate from the two Feynman 
diagrams in subprocess 2.2, W-gluon fusion: tqb . These calculations were done in 

the unitarity gauge with the CTEQ3M parton distributions. The Feynman diagram in which 
the initial state gluon produces a bb pair contributes 145% to the total rate, whereas the 
diagram in which the gluon splits to a it pair contributes only 5%. Although this diagram 
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contributes so little to the total, the destructive interference between the two diagrams is 
very significant, at -50%. Figure 3 shows the contributions to the total W-gluon fusion 
cross section from each diagram and from the interference between them as a function of 
the top quark mass. 



4. Wtb Coupling and V tb 

We have investigated the effect of an anomalous contribution to the Wtb coupling, in 
the form of a right-handed (V+A) structure, specified by a parameter A r . In the unitarity 
gauge, the Wtb coupling is given by: 

where e is the electronic charge, sim} w = 0.474 and y '„ and y 5 are Dirac matrices. 

The dependence of the total single top cross section on the parameter A r is shown in 
Fig. 4 for the upgraded Tevatron energy of s[s = 2.0 TeV and V tb = 0.999. Here 
o(t + i) = a((tb) + (tb) + (tq) + (iq) + (tW~) + (tW + )) . The cross section rises from the SM value 
of 2.37 pb to 4.60 pb when A r = +1. 
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Fig.4 Total single top quark production cross 

section in pp interactions as a function of 
the right-handed (V+A) coupling parameter 




0.5 | .... | .... | .... | .... | 
-1 -0.5 0.5 1 

(V+A)/{V-A) Coupling Strength A 

Fig. 5 Region of sensitivity at the 95% C.L. in the 
( V,b,A r ) plane, for an experiment with 2 ftf 1 

A,. of data, S:B = 1:2, 10% efficiency and m, = 
180 GeV. 



We have estimated the experimental sensitivity in the (V tb ,A r ) plane from the next 
Tevatron collider run. The parameters used in our estimate are that the collider will operate 
at 2.0 TeV, there will be a 2 fb 1 data set, the top quark has a mass of 180 GeV, and the 
search will take place in the electron+jets and muon+jets decay modes with 10% overall 
efficiency. We assume the branching fraction will remain constant at 2/9 and that the 
signal-to-background ratio will be 1:2, based on a study for the TeV-2000 project 14 . If the 
Wtb coupling has SM form and there are only three quark families (i.e. 
0.9988 < V A < 0.9995 at the 90% CL 15 ), 
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then there should be approximately 105 fully reconstructed b tagged events. Assuming that 
315 events are seen (105 signal and 210 background), limits can be placed on V tb as a 
function of any anomalous coupling as parametrized by A r . These 95% CL limits are 
shown in Fig. 5. If the W tb coupling is purely left-handed, then the limits on V lh are 
0.817 <V tb <l. For A r = 0.5, 0.735 <V tb <\, and if A, = 1.0, then 0.586 <V tb < 0.826 . 
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